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rZ. t The variational principle of barotropic Eulerian fluid dynamics is known to 

O ' be quite cumbersome containing as much as eleven independent functions. This 

*C/2 ! is much more than the the four functions (density and velocity) appearing in 

^"v the Eulerian equations of motion. This fact may have discouraged applications 

'-~h | of the variational method. In this paper a four function Eulerian variational 

principle is suggested and the implications are discussed briefly. 

1 Introduction 

O 

The motion of a fluid is usually described in two ways. In the Lagrangian approach 
we follow the trajectory of each fluid particle. While in the Eulerian approach we 
^ ! study the evolution of the velocity and density fields. 

q^ | Variational principles have been developed in both cases. A variational principle 

in terms of Lagrangian variables has been described by Eckart [y and Bretherton 
[ |] . While a variational principle in term of the Eulerian variables is given by Herivel 
[|, Serrin @j, Lin § and Seliger & Whitham f|. 

The Eulerian approach appears to be much more appealing involving measurable 
Q_|. quantities such as velocity and density fields instead of the trajectories of unseen 

particles of the Lagrangian approach. 

However, regretfully the variational principle of the Eulerian flow appears to be 
much more cumbersome than the Lagrangian case, containing quite a few "Lagrange 
multipliers" and "potentials". In fact the total number of independent functions in 
this formulation according to the approach suggested by Herivel |J, Serrin Q and 
Lin H is eleven which exceeds by many the four functions of velocity and density 
appearing in the Eulerian equations of a barotropic flow. 

The variational principle of the Lagrangian approach is on the other hand simple 
and straightforward to implement. Bretherton [2[ has suggested to use the La- 
grangian variational principle with Eulerian variables and especially constructed 



"Eulerian displacements" for which the equations of motion are derived. This pro- 
cedure appears to be inconvenient for stability and numerical calculations. 

Instead I intend to develop in this paper a simpler variational principle in terms 
of Eulerian variables. This will be done by rearranging the terms in the original 
variational principle and using a simple trick to be discussed below. The result will 
be a variational principle in terms of four Eulerian independent functions which is 
the appropriate number. 

The plan of this paper is as follows: first I review the original Eulerian variational 
principle and the equations derived from it. Than I derive a new variational principle 
and discuss its implications. 

2 The Eulerian Approach to the Barotropic Flow 

In the Eulerian description of barotropic fluid dynamics we seek a solution for the 
velocity v = v(x k ,t) and density p = p(x k ,t) fields of a given flow. Those fields 
depend on the space coordinates r = (x k ) = (x, y, z) [k, l,m,n = 1, 2, 3] and on the 
time variable t. The fields are obtained by solving the Euler equations: 

dv -> -> 

— + v-W+V(/i + $) = (1) 

in which V = -^. The potential <I> = &(x k ,t) is a given function of coordinates, 
while h = h(p) is the specific enthalpy which is a given function of p usually des- 
ignated as the equation of state. In addition the density and velocity fields must 
satisfy the continuity equation : 



g + v.« 



0. 



(2) 



Thus we have to solve four equations in order to obtain four unknown fields. This 
can be done when supplemented the appropriate boundary conditions. 



2.1 The Classical Eulerian Variational Principle 

Researchers (Herivel H , Serrin H , Lin || ) seeking a variational principle from which 
the Euler and continuity equations (equation (d) and equation (g)) can be derived, 
arrived at the action A given in terms of the Lagrangian L as: 

rtx 

A = Ldt with 



L 



to 



V 



V 



1^ 



e{p) - $ 



pd 3 x 



,,| + v.(h)-^-^ 



d A x. 



(3) 



In which s(p) is the specific internal energy connected to the specific enthalpy by 
the equation : 

d(pe) 



h 



Op 



(4) 



The volume element is given by d?x and the operator -^ is defined such that: j^ = 
§l+v-V . The variation principle contains in addition to the desired four functions p, v 
also the seven "potentials" is, a, (3 reaching to a total of eleven variational variables. 
Varying the above action with respect to those variables such that the variations 
vanish for r and t sufficiently large. We see that in order to have 5A vanish for 
otherwise arbitrary (5v,dp,5v,8a,5/3), the following equations must be satisfied: 

^ + V-(pO)=0 (5) 

^ + ?-G**)-0 (6) 

$ = a-V/3 + Vv (7) 

f - • 

Combining equation @ and equation (JSj) we arrive at the simpler equation : 

^ = 0. (10) 

Dt v ; 

Calculating the expression j¥ using equations ^gjl^ii) we arrive at the equations 

of Euler: 

Dv dv -» -> 

— = - + v-Vv = -V(h + ^. (11) 

Thus the functions v and p are extrema of the action A if they satisfy the Euler and 
continuity equations . 

2.2 The Variational Principle of Seliger &z Witham 

Seliger & Witham || have proposed to take a. = (a, 0,0) and j3 = (/?, 0, 0), in 
this way one obtains a variational principle with only seven functions. Rewriting 
equations (|5[0,^J9,1O) for this case we obtain the following set of equations : 



% + v-(H = o (12) 

V = aV/3 + W (13) 

^ = (14) 

Dt y ' 

D£ 
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Dv 1 



(15) 

v 2 - h - $. (16) 



Dt 2 

We see that equations (|12|, jig ) remain unchanged. The velocity v is given now by the 
Clebsch representation (equation ( |l3| ) ) and the entire set of equations are designated 



as Clebsch's transformed equations of hydrodynamics (Eckart ffl, Lamb m p. 248). 
Bretherton H have quoted a remark by H. K. Moffat concerning the possibility of 
the function v being not single valued. According to Moffat || the helicity integral: 



H 



V x v ■ v d x 



(17) 



which is a measure of the knottedness of the vortex lines, must be zero if v is single 
valued which is clearly not true in general. This can be shown as follows: suppose 
that v is single valued than by inserting equation (|l3|) into equation (|l7|) and 
integrating by parts we obtain: 



n 



Vi/-Vax Vf3d 3 x = / V • (i/Va x Vf3)d 3 x = / iA7 x v ■ dS 



(18) 



If we choose a volume made of closed vortex filaments such that Vxf dS = than 
clearly TL = 0. However, as Bretherton Q noticed the equations of motion being 
local are unaffected by global properties such as the non single- veluedness of v. 

An analogue from classical mechanics may make things even clearer. A few 
will object to describe the two dimensional motion of a particle moving under the 
influence of the potential V(R) (where R is the radial coordinate), by the coordinates 
R, (p. Where (j) is the azimuthal angel which is not single valued. This form is found 
to be more convenient than the single valued Cartesian x, y representation. 



3 The Reduced Variational Principle 

Although Seliger & Witham (1968) have managed to reduce the Eulerian variational 
principle from eleven to seven functions. The number of variational variables is 
still too much, since the Eulerian equations contain only four unknown functions. I 
intend to suggest a solution to this problem here. 

First let us rewrite the Lagrangian of Seliger & Witham (equation (y) with a 
and j3 being scalars): 



L = 

+ 


I 


-v 2 - e{p) - $ 


pd 3 x 
)^-pa(^ + v-V(3) 


Next we rearrange terms: 


L -i 


d/3 dp' 

- pa Tt + v Yt 


d 3 x 


+ x 


'1 ^2 

2 PV 


+ uV ■ (pi 


1)-pc 


w- V/3 


d 3 x - f [e(p) 
Jv 



d s x. 



(19) 



(20) 



Furthermore, we introduce the identities: 



dp 

uV ■ (pv) 



d(vp) dv 
~dt P ~dt 
V • (vpv) — pv • Vv. 



(21) 



Inserting the above identities into equation (EQ) and rearranging terms again we 
have: 



L 



dp dv 



pd 3 x 



-pv 



pv ■ (aV/3 + Vv) 



d 6 x 



[e(p) + $] pd 3 x 



^d 3 x + J v-wpruf:,: 



Now since: 



]-v 2 -v-(aVP + Vv) = hv - aV P -Vvf -haV p + Vvf 



We finally obtain: 



in which: 
And 



, 8(5 dv, 



L v 



L — L r + L/ v 



(aV/3 + Vvf - e(p) - $ 



pd 3 x + 



dt 



-{aVP + Vu-vypd 3 x+ / V • {vpv)d 6 x. 
v 2 Jv 



(22) 

(23) 
(24) 

(25) 
(26) 



The Lagrangian L is dissected into a part depending on the fluid velocity v which is 

denoted as L v and the remaining part which does not depend on v which is denoted 

L r . Let us look at the action: 

r ti 

A, r = / L r dt. (27) 

J to 

This action is a functional of the four variables: a, (3, u, p. Note that the last term of 
L r is a full time differential and thus does not contribute to the equations of motion. 
In order to simplify the variational calculations we introduce the notations: 



-n - , D d _ -* 

u = aV 3 + \v and — — = — + u ■ V. 

Dt at 



(28) 



Using the above notations and assuming that the arbitrary variations (5p, Sis, 5a, 53) 
vanish for f and t sufficiently large we obtain that 5A r = only if the following 
equations are satisfied: 



| + v-« = o 



d(pa) 

~dT~ 

Dp 
~D~i 
Dv 
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+ V • {(ma) = 





1 



h-$. 



(29) 
(30) 
(31) 
(32) 



CombiningF] equation ( |29] ) and equation (BO) we arrive at the following set of 
equations : 

^ + V-(/m) = (33) 

^ = (34) 

Dt y ' 

# = (35) 

Dt K ' 

El = U 2 -h-<S>. (36) 

Dt 2 y ' 

If we take u to be a Clebsch representation of some velocity field v = u, than the 
above equations become the Clebsch transformed equations of fluid motion (equa- 
tions ^2 . 14,15,16J). Thus we have achieved a variational principle that does not 



contain v as a variational variable. This situation in which only the potentials ap- 
pear in the variational principle but not the physical velocity field v itself is familiar 
from other branches of physics and will be discussed below. 

It remains to discuss the Lagrangian L v given in equation (£6|) , using the theorem 
of Gauss and the definition (^) of u we arrive at the form: 

L v = f ~(it- v) 2 pd 3 x + J vpv- dS. (37) 

If we take u to be a Clebsch representation of v than the first part of L v vanishes. 
Further more the surface integral remaining vanish for both fluids contained in a 
vessel in which v is parallel to the vessels surface and for astrophysical flows in 
which the density vanishes on the surface of the flow. However, in the case of helical 
flows v is not single-valued and one must add in addition to the physical surface a 
"cut". Thus one obtains: 

L v = J [v]pv-dS (38) 

J cut 

where [i>] represents the discontinuity of v across the "cut". We conclude that 
for helical flows L v does not vanish altogether. And thus for those flows the new 
Lagrangian L r is not equal to the Seliger & Witham Lagrangian L, the difference 
being a surface integral over the "cut". 

3.1 A Comment Regarding Esthetics 

It was already noted that the last term of L r given in equation (^5|) does not 
contribute to equations of motion. It seem esthetically plausible to omit this term 
altogether thus obtaining: 



L. 

iv 



{a^ + ^)-\{aV(i + Vvf-e(p)-^ 



pd 3 x. (39) 



The Lagrangian is quadratic in the spatial Clebsch representation aV/3 + Vz^, and 
linear in the temporal "Clebsch representation": a-M- + — 



dt ~ dt- 



1 Compare the transition from equation (pi) to equation (|l0|) . 



3.2 Potentials as Lagrangian Variables 



The Lagrangian L s contains variables (a, /3, v) that have no physical meaning outside 
the context of the Clebsch representation. This is also the situation in electromag- 
netics in which the Lagrangian (see for example Goldstein (1980) p. 582) is given 
by: 



L, 



v 



L.dA 

2"~di 



^(^T +2VA -^ + VA 2 



dA 
~dt 



\2 



V x A 2 ) + e(j ■ A - pA c 



d 3 x. (40) 



In which A, Aq are the vector and scalar electromagnetic potentials, e is the charge 
of the electron, j is the current and p is the charge densityf]. Only after the problem 
is solved in terms of the potentials A, Aq the physical electric E and magnetic B 
fields can be obtained through the relations: 



E 



9A A A 

■m~ VAo 



B = V x A 



(41) 



which are analogue to the Clebsch representation given by equation (O) in the fluid 
dynamical case. 

3.3 Linearity in Time Derivatives 

Another odd characteristic of the Lagrangian L s given by equation ( |39|) is that 
it contains only linear terms in time derivatives. This is unlike the generic case 
in classical mechanics in which terms quadratic in time derivatives appear in the 
kinetical energy part the Lagrangian. For this situation I have been able to find an 
analogue from quantum mechanics. The Schroedinger Lagrangian is given by: 



L 



sell 



Bib h 2 - 

-i?i^-ip* + W • V^* + Vipip* 

at 2m 



(Tx. 



(42) 



in which h = ■£- is the Planck Constant, % = \J — 1, ifi is the wave function and tp* is 
its complex conjugate. L sc h is linear in time derivatives and in this respect resembles 



4 Conclusions 



Although a compact variational principle have been obtained in equation ( p9| ) it 
remains to see the possible implications of this expression. One possible utility is 
to obtain better numerical algorithms for solving fluid flow problems based on the 
above variational principle (see Yahalom [|10[ ) . Another may be to study the stability 
of certain stationary flows in the spirit of the works by Katz, Inagaki & Yahalom 



11] and Yahalom, Katz & Inagaki [12|. 



J Not the matter density! 



Further More, it is also desirable to have a variational principle for incompressible 
flows which should contain only three functions since the density is not a variable 
in this case. A similar reduction in degrees of freedom should be obtained for two 
dimensional flows in which the velocity has only two components. The above list 
of extensions of the Eulerian variational principles is of course not comprehensive. 
And the undertakings mentioned are left for future papers. 
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